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In memory of Leonard Carlitz 
Abstract 

The exponential generating functions ^ n>1 n n_m ^7 (for m 6 Z) are proven 
to be rational functions of the tree function CP, ^. n" -1 ^) and then of the end- 

ofunction function (X}»>1 ""^r)) an< i some coefficients of these rational functions 
are identified with certain generalized Stirling, Eulerian, and Bernoulli numbers. 
Of the four families we cover, one had been previously analyzed by L. Carlitz. 



1 Introduction. 

In Jy, Carlitz expressed J2 n >i n ™ _m fn" ^ or ^ — 171 e ^ as a rational function of A(= 
^( z )) — En>i n"~ 1;5 T- The coefficients in this rational function (written in lowest terms 
with factored parts) are positive integers which had previously appeared in H and were 
subsequently given combinatorial interpretations in 1 1 and These numbers are now 
known as "second-order Eulerian numbers" (cf . M ) . 

In this note, we consider to > 1 and find polynomials in A. We show that the 
coefficients of these polynomials are the virtual Stirling Numbers of the First Kind 
studied by D. Loeb in Q. Previous tabulations and interpretations had appeared in Q 
and elsewhere in the "finite differences" literature. 

We introduce the variable ( — X)n>i n "frr = T^T an< ^ investigate the coefficients 
in the rational functions of £ equivalent to those in the above sequence. In the case 
TO < we find polynomials of degree 2m + 1 which, in factored form, display integer 
coefficients equal to the "associated Stirling Numbers of the second kind" studied by 
Riordan in pj§ . These have also been named "second-order Stirling Numbers" (in ||, 
for instance) and we show that easy manipulations produce new identities relating them 
to the second-order Eulerian numbers. 

It remains to address the £ expressions in the case to > 1. The rational functions, in 
lowest terms, are quotients of two polynomials each of degree |m| in £. While we obtain 



satisfactory results about the asymptotic behavior of the coefficients, the question of 
combinatorial interpretation is open. 

2 An Integral Operator. 

For integer values of m we define the sequence G m (A) := R m (z) := J2 n >i n " _m frr- 
Thus, as defined above, A = Gi(A) and ( = Go (A). We define the formal integral 
operator I by I(F(A)) = F m (p)±=* dp. 

Theorem 1. I(G m (A)) = G m+ i(A). 
Proof. 

pz du f 1 p f 1 p 

G m+1 (X) = Rm+^z) = / R m (u) — = / R m {pe~ p ) dp = / G m (p) dp □ 

Jo u J Q P Jo P 

Corollary 1. G m _i(A) = j^G'JX). 



Remark 1. The differential operator in the Corollary is essentially that used by Carlitz 
Specifically he derived the following transform for nonnegative m: If 

n . m 

n>l y ' k=0 

then 

Here the symbol (( ™ )} has become known as the second order Eulerian number Ta- 
ble 270]. Combinatorial interpretations o/((^)) were subsequently given by Riordanpl\] 
and Gessel and Stanley j^j. 

Corollary 2. If m > 0, G m (A) is a polynomial of degree m without constant term 
(Appendix 1). If we write G m (A) = X)fc=i ( — l) fe ~ 1 f ( m > k)^ k r then, form > 1, 

k ■ g(m, k) = g(m — 1, k) + g{m — l,k — 1). (3) 

We note that g(m, 1) ~ 1 and set g(m, k) = for k > m > 0. These two boundary 
conditions, together with (3), determine g(m,k) for m, k > 1. A combinatorial inter- 
pretation of g(m, k) is as follows: let N = m — k + 1 and consider the set of all k x N 
matrices such that each column is a permutation of {1, 2, . . . , k}. Then g(m, k) is the 
probability that, in a randomly selected matrix from this set, no one number is below 
the number 1 in every column. 
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We note two previous appearances of the numbers g(m, k) in the literature. In [g, 
Appendix] , a definition of generalized Bernoulli numbers of negative degree is given in 
accordance with which g(m, k) — — ( m ^ 1 ^j BL^' . In the Stirling numbers of the First 
Kind are generalized such that g(m, k) = (— l) fe s(— k, m — k). Many formulae for g(m, k) 
may be found in these references and elsewhere (cf. Appendix 3). 

3 Pullback by Fractional Linear Transform (m < 0) 



If A may be called the tree variable, then 

n>l 

is the endofunction variable. When we express R m { z ) = G m (X) in terms of C a beautiful 
phenomenon occurs: Carlitz' rational functions (m nonpositive) become polynomials in £ 
with integer coefficients strongly related to numbers called "associated Stirling numbers 
of the second kind" by Riordan [R2] or "second-order Stirling numbers" (denoted {{ ™ }}) 
by Fekete[F]. We adopt the latter notation, in terms of which {{ ™ } } = k ■ {{ m ^ k 1 }} + 
(m — 1) • {{™~i 2 }} is the defining 3-term recurrence relation. 



Proposition 1. For m G N. 

G- m (X) = (l + Cr +1 E{{ m + fc }}c fe :=H_ m (0- (4) 

Proof. Making the substitution A = in (1), and writing the Carlitz differential 
operator as £(1 + C) 2 ^> w e calculate, with G_ r „(A) = (1 + () m+1 Y^k=i a ( m > k)( k and 
Theorem 1, that a(m, k) = (m + k— 1) • a(m — 1, k — 1) + k ■ a(m — 1, k), from which the 
result follows. □ 

These formal calculations have some pleasant concrete consequences. 
Corollary 3. 

£<G><i + <r-*v-£;{{T}}<" < 5 » 

fe=0 fc=l ' ' 

Summation identities are produced by specializing ^. A pure counting identity follows 
from the binomial theorem: 

Corollary 4. 
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This is completed to an inverse pair of identities (as in p^|) by 
Corollary 5. 

PHB-irt;:; 1 )!!":;: 1 }} < 7 > 

4 Pullback by Fractional Linear Transform (m > 0) 



The C~ expressions for positive m may be written H m {Q = (1+C)~ m Sfeli M m ; k)( k . 
The rational numbers /i(m, fc) do not appear, to our knowledge, in the literature (see 
Appendix 2, top). In particular, combinatorial interpretation of h(m, fc) is unknown. It 
is possible, of course, to use our operators and the correlative properties of g(m, fc) to 
analyze these coefficients quantitatively. We state the basic facts without proof. 

First we note that these coefficients may be given by the formula 

h(m,k) = p^(-iy- 1 (^2-)9(m,j) (8) 

which may be combined with known formulas for g{m, fc). A better clue toward finding 
a combinatorial interpretation may come from the rather simple recurrence relation 

fc • h(m, fc) = h(m - 1, fc) + (m - k + 1) • h(m, fc - 1). (9) 

Since g(m, fc) has an easy interpretation as a probability, one might hope to interpret 
h{m, fc) as an expected value in some finite problem. To this end we note that both 
l!2! • • • ml ■ h(m, fc) and (l!2! • • • fc!) m ~ fe+1 • h(m, fc) are integral. We further state 

m 

J2(-iy- 1 Krn,j) = l/m\ (10) 
i=i 

These facts may be easily proved by induction. It is also easy to verify that 

lim h(m,m) — l/e (11) 

m — >oc 

The growth of /i(m, fc) as m — > oo for fixed fc is of interest. Experiment strongly 
suggests that (fc — 1)! • h(m, fc) is asymptotic to a polynomial pk{m) of degree fc — 1 in 
to with nonnegative integer coefficients. 

5 Conclusion. 



The variables A = T(z) and £, as well as W := — T{— z) have been used extensively in 
analyses of random mappings and random graphs (se j?], |4j for surveys) in which other 
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R m {z) appear. The latter also occur naturally as exponential generating functions of 
tree-like structures. We hope that our very basic results may be useful to further such 
studies. 
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Appendix 1: The A Sequence 

For A(= T{z)) = X)„>i n " _1 fn" an d integer values of m we define the sequence 
G m {\) := R m (z) := £ n n - m ^'. Then 



Ge{X) - "720 A + 7200 A ~ 3456 A + 1296 A ~ 32~ A +A 



g ^= ik A5 -S A4+ ^ A3 -i A2+A 



g 4 (a) = -1a 4 + 11 A 3_ 7 a2 + a 



G 3 (A)= 1 A 3_ 3 A2 + A 

G 2 (A) = -^ A2 + A 



Gi(A) = 
G (A) = 

G-i(A) = 

G_ 2 (A) = 



G_ 3 (A) = 
G_ 4 (A) = 
G_ 5 (A) = 



1-A 

A 

(1~AF 

A(l + 2A) 
(1-A) 5 

A(1 + 8A + 6A 2 ) 

A(1 + 22A + 58A 2 + 24A 3 ) 
O^Ap 

A(l + 52A + 328A 2 + 444A 3 + 120A 4 ) 

0~W 
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Appendix 2: The £ Sequence 

For £ — X)„>i n "frr an d integer values of m we define the sequence H m (() 
Rm(z) :=E„>i« n - m £- Then 



#5«) 



#4«) 



^3(0 = 



MO 



/7l -I- 129 A 4- 8513 r 2 4- 691 A3 I 96547>4 , 96547 a5\ 
"•A 1 32 ' 1296 ^ ^128^ ' 43200 ^ ' 259200*' > 

(i + C) 6 



f(\ 4- 19/ 1 4- 1547 a2 , 1631 a3 , 1631 M\ 
Si 1 T 16 <, 4- 432 s 4- 864 <■, -r 43 20 ^ ) 

(i + C) 5 



c(i + ±fc + fc 2 + ^c 3 ) 



(i + C) 4 

c(i + |c + Ac 2 ) 



(i + C) 3 

co + io 

(i + C) 2 

c 



i + C 

#o(C) = c 

ff-i(0= C(i + C) 2 

ff- 2 (0= c(i + C) 3 (i + 30 

fl-_ 3 (0= C(i + C) 4 (i + 10C + 15C 2 ) 

H-i(C) = C(l + C) 5 (l + 25C + 105C 2 + 105C 3 ) 

ff_ 6 (C) = C(l + C) 6 (l + 56C + 490C 2 + 1260C 3 + 945C 4 ) 
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Appendix 3: Extended Remark on g(m, k) 



The tabulation of ™ led Carlitz to ask for combinatorial interpretations. In con- 
trast, the coefficients g(m, k) are well known. In 0, a definition of generalized Bernoulli 

numbers of negative degree is given from which we see that g(m,k) = -(^^B-r 
Also g(m, k) — (— l) fe s(— fc, m — k), a virtual Stirling number as defined in 

Consider the recurrence relation 

k ■ g(m, k) — g(m — 1, k) + g(m — 1, k — 1) 

for m, k > 1 where we initialize by g(m, k) — if k > m, and <?(m, 1) = 1. Since this is 
our definition, the verification of the following formulas is routine. The unique solution 
may be expressed as a sum of Egyptian Fractions: 

g( m> k) = - ' 



l<il <%2< - ■ -^m — fc</C 



An alternative form is related to an iterated difference of f(x) — xm lk+i '■ 

. fc-i 



'' 1 f k - l\ (-]}> 



,._„ v J J{j + l) m - k+1 

From the recurrence for the g(m, fc)'s also follow generating functions: 



g(m,k) = [z m - k ] f[ — 



v-z 
P =i 

which may be written in a connection coefficient formula: 

1 



^2(-l) n g(k + n,k)x n 



(x + l)(x + 2)---(x + k) n 
(linking g(m,k) with the signless Stirling cycle numbers). 

An integral of a polynomial over the hypercube yields the same numbers: 

g(m, k) = 1 / [ ■■■ [ (1- X!X 2 ■ ■ ■ x m - k +i) k ^ dxi dx 2 ■ ■ ■ dx m ^ k +i 
{k - 1)! Jo Jo Jo 

since it evaluates to the iterated difference form above after applying the binomial the- 
orem. 



The rational number g(m,k) has the natural denominator (kl) m k+1 . The corre- 
sponding numerators N(m, k) form a triangle of integers generated by 



N(m,k) = [z m - k ] 



-z 



and previously tabulated and studied in David, Kendall, and Barton, Symmetric Func- 
tion and Allied Tables, Cambridge Univ. Press, London, 1966 (Table 5.4.1). 
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Web Resource 

Some additional data on the questions raised in Section 4 may be found at 
http : / / www . saturn . math . uaa . alaska . edu/~smiley/numbertriangle . html 
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